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The time independent mean field method (TIMF) for the calculation of matrix elements of non 
relativistic propagators is based on a variational principle whose non linearity induces a multiplicity 
of variational solutions. Several of them can break any symmetry shared by the Hamiltonian and 
initial and final states. We describe a soluble model where, in particular, time reversal and parity 
breakings occur. Such breakings account for important properties of propagation amplitudes. 
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The observation of symmetry breakings in solutions 
of variational principles is familiar for static problems, 
such as those described by the Hartree-Fock method 
for deformed nuclei. Similar breakings can be expected 
for other variational theories, like the time independent 
mean field (TIMF) theory of reactions and decay . The 
TIMF method has already been implemented for various 
nuclear and atomic processes involving 3, 4, 5 and 8 parti- 
cles and symmetry breakings were indeed found |^ . In 
\ this letter we describe a soluble model which completely 
' interpolates between mean field trial states and exact so- 
, lutions. It provides a detailed observation of symmetry 
■ breakings as well as conservation. Because propagation 
amplitudes are obtained by saddle points rather than ab- 
solute maxima or minima, the multiplicity of solutions 
raises an ambiguity. The model shows that actually all 
' ■ mean field solutions should later be mixed linearly to re- 
O , cover, for the complete energy range, the best possible 
5 ■ approximation of the exact solution. 

We start from a Schwinger-like functional [|1|,| 
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to calculate the amplitude D{E) =< x'\{E - H)~'-\x > 
for propagation between initial and final states XjX' ■ The 
variation of F gives, with appropriate normalizations, 

{E-H)\^S/> =|x>, (2a) 
<^'\{E-H) = <x'\, (2b) 
DiE)^<x'\^> =<*'|x>, (2c) 

with (^f, ^P') as a saddle point of F. For our model we con- 
sider a single-particle basis made of two square integrable 
orbitals \a > and \b >, with single-particle energies —1 
and 1, respectively. In the resulting space for two distinct 
particles, spanned by |a^ >, \ab >, \ba >, >, we use for 
H simply the corresponding sum of one-body Hamiltoni- 
ans, with eigenvalues —2, 0, and 2, respectively. This is 
clearly a Hamiltonian invariant under particle exchange 
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(denoted parity in the following) and complex conjuga- 
tion (time reversal). Then we choose |x' >= \x >> and 
X as a two-body state, described by a product of two 
identical, single-particle mixtures, 

|X>= (r|a>+|6>)i(r|a>+|6>)2, (3) 

where r is a frozen, real parameter. Both x, x' are invari- 
ant under parity. The same holds for ^ . Furthermore, 
since |x >= |x' > and r is real, time reversal symmetry 
is expressed by the property >= 1^* > if is real, 
and l^f' l^** > if is complex, sec Eqs.(2a,b). 
It is seen trivially that the exact value of D{E) is 

^4 1 
D-(^) = ^, + ^ + ^,- (4) 

This pole structure allows E to remain real in the fol- 
lowing. Our model investigates reconstructions of D{E) 
with trial functions ^, ^' restricted to less flexible forms, 

1$ > = c| > -I- s I aib2 > + t\ bia2 > + st \ b'^ >, 
<^'\ = c<a'^\+u< aib2\ +v < 6102] +uv < (5) 

where c is a frozen, real parameter which expresses corre- 
lation constraints, and s, t, u, and v are four variational 
parameters, possibly complex. It is easily verified that, 
for the special value Cgx = E"^/ {E^ — 4), the space of trial 
functions $, contains the exact solutions 4", 4"' given 
by Eqs.(2), with s = t = u = v = E/[r{E - 2)]. Con- 
versely, for c = 1, we generate a TIMF limit where the 
trial states $, arc uncorrclatcd, namely simple prod- 
ucts of single-particle orbitals. The model thus achieves, 
as a function of c, an interpolation between a mean field 
approximation and an exact theory. Note that, despite 
its simplicity, the model is non trivial since the inverse 
of an additive one-body operator is a complicated sum of 
many-body operators. This letter studies whether time 
reversal, namely s = u and t = v, on one hand, and on 
the other hand parity, namely s = t and u = v, are prop- 
erties conserved or broken by the variational principle. 
With the restricted ^, $' of Eqs.(5), F becomes 

p {cr^ + rs + rt + st){cr^ + ru + rv + uv) 

~ E{c^ + su + tv + stuv) - 2{stuv - c^) ' ^' 

whose derivatives with respect to s, u, v are easy to ob- 
tain. Let Fg, Ft, Fu, Fy be the multivariate polynomial 
numerators of these derivatives. From the conditions 
Ft = Fy = Oa, straightforward, linear elimination of t and 
V in terms of rational functions of s and u is possible. The 
substitution of such rational fractions for t and v inside 
Fg and Fu then leaves two residual polynomial conditions 
for s and u. It is finally easy, albeit slightly tedious, to 
eliminate for instance u between these two polynomials 
and obtain a polynomial resultant R{s; E, r, c), the roots 
of which designate the saddle points of F, Eq.(6). 
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It turns out that, except for spurious factors which 
either identically cancel F or make divergent, R 

reduces into the product of three simpler resolvents, 

Rcc = c^r{E + 2) + c(2c + cE - Er'^)s + {E + 2cr'^ - cEr'^)s^ + r(2 - E)s^, (7a) 
Rpb = cE{2c + cE + Er^) + r{Ac^ + E'^ - c^E^)s + E{E - 2cr^ + cEr^)s^, (7b) 
Rtb = cr{2c + cE + Er^) + c(2r^ - Er'^ + E + 2)s + r{2cr^ - cEr^ - E)s^. (7c) 

The first factor, R^c, generates four saddle points for 
which s = t = u — V, namely both parity and time 
reversal are conserved. The second factor, Rp^, induces 
two saddle points for which s = u and t = v, but s ^ t, 
hence parity is broken, while time reversal is conserved. 
The corresponding value Fpb of F is the same at both 
such saddle points, 

[r^ic^r^ + c2 + 2c - 1) + c^]E^ + 2c^l - r^)E - Ac^r'^ 



E{c'^E^ _ £;2 _ 4c2) 

(8) 

The last factor, Rtb, induces two saddle points for which 
s = t and u = V, but s ^ u, hence parity conservation 
but time reversal breaking, with again equal values of F 
at both saddle points, 

{cr^ + 2r^+c)E + 2cil-r^) 
Ftb = {c-l) (c2_ 1)^2 _ 4,2 • (9) 

While many combinations of the real parameters E,r,c 
induce real values for s,t,u,v at such symmetry break- 
ing saddle points, see Figs. (1,2) for a pb- and a tb saddle 
point, respectively, there are equally frequent situations 
where those pairs (s,t) or (s,m) which are not degener- 
ate actually take complex conjugate values. Neverthe- 
less, as shown by both Eqs.(8,9), the functional remains 
real valued. This allows useful analytic continuations of 
variational estimates of _D, which will indeed be plotted 
in several forthcoming Figures, regardless of the real or 
complex nature of the corresponding saddle points. Con- 
versely, for the fully symmetric cases described by Rcc, 
real estimates of D are obtained when s is real only, un- 
less ReF is used as an estimate of D when F is complex. 

By means of Fig. 3 we now discuss the behaviors of 
the six estimates of D provided by R when c takes on 
all values between the limit without correlations, c — 1, 
and the limit allowing an exact result, c = Cex- We set 
r = —1.2 and E = 1.5 for Fig.3. Hence Cex = —1.29 and 
Dex = 0.51. The full lines show two branches of estimates 
provided by Rcc- Its other two branches turn out to be 
either complex or out of range in that case. Not surpris- 
ingly, one of the full lines contains the point {cex, Dex), 
representing the fully symmetric, exact solution ^, 

It is remarkable that the long-dashed line, which cor- 
responds to the roots of Rtb, also contains this exact 
point. It is also remarkable that the time reversal break- 
ing, measured by |s — w| for instance, stays finite on that 
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long-dashcd branch when c approaches c^x and does not 
vanish at that point. It is nevertheless easy to verify that, 
when c = Cex, then Ftb = D^x as a function of E and r. 
This indicates that, as long as correlations are allowed 
in the trial states, even partly only, the breaking of time 
reversal symmetry may help in finding good estimates of 
matrix elements of propagators. The paradox of a finite 
value of s — It at the exact point is understandable, be- 
cause a well-known property of F, see Eq.(l), is that only 
any one of the two equations, Eqs.(2), needs to be solved 
exactly to obtain Dg^- It is unfortunate that, as shown by 
Eq.(9), the stationary value Ffb provided by uncorrelated 
trial states, namely for c = 1, identically vanishes. 

The second full line branch goes through that same 
point {c,D) = (1,0). This is no accident, because i?cc 
gets a trivial, energy independent root s = ~r when c = 
1, for which the fully symmetric functional identically 
vanishes. Only three fully symmetric saddle points are 
left to depend on E for uncorrelated trial states, together 
with the pb branch. In the present case, besides that 
trivial fourth fully symmetric root s = — r = 1.2, one fully 
symmetric root is real and the other two are complex. 

The short-dashed line on Fig. 3 represents the behavior 
of the pb branch. It shows a vague trend to interpolate 
between the two real, fully symmetric branches, but its 
utility will be apparent later only. Indeed, it does not go 
through the exact point {Cex,Dex), nor does it show an 
accurate reproduction of D^x for c = 1. 

We show on Fig. 4 the multivalued nature of the esti- 
mate(s) provided by the four roots of Rcc- Obviously, the 
number of real roots depends on E, r, c. 

We now turn to behaviors of the estimates of D as func- 
tions of E. This is illustrated by Fig. 5, where Dex{E) (full 
line) is compared to its various estimates with c = 1. It is 
seen at once that a fully symmetric estimate (dashed line) 
is convenient to reproduce pole behaviors near E = ±2. 
This is not completely siuprising since the correspond- 
ing eigenstatcs \a? > and \b^ > are themselves invariant 
under particle exchange. Conversely, eigenstates such as 
ai&2 > a-nd |6iCi2 > are not invariant under such an ex- 
change, and it is thus reasonable that a pb estimate is 
needed to account for the pole at E = 0. Moreover, it is 
consistent that the pb estimate fails for E ~ ±2, and the 
fully symmetric estimate in turn fails for E ^ 0. For the 
sake of completeness, we also show on Fig. 5 the behavior 
of the tb estimate obtained when c = 0. It is clear that 
this tb estimate, like the pb one, succeeds for ~ and 
fails for E ~ ±2. All these properties are transparent 
from the corresponding special values of the functionals, 




{Er"^ + E + 2)(2r2 - Er"^ - E) 
4B 



(10a) 




2r^ 



(10b) 





(10c) 



{s'^ + l){Es'^-2s'^ + E + 2y 
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The residues of the special values, Eqs.(10a,b), of Fpf, and 
Ftb at E = are and 2r^, respectively. This compares 
well with the exact residue 2r^, see Eq.(4). The three 
non trivial saddle points of Fcc|c=i are the roots of 

i?cc|c=i = r{2 - E)s^ + Es'^ - Ers + (2 + E). (11) 

A straightforward elimination of s between Eq.(ll) and 

Eq.(lOc) generates a cubic polynomial condition relating 
these three saddle point estimates Dec to r and E, 

16(4 - E'^)Dl^ + 8[12 + QE- 2E^ - E^ + 6(^2 - A)Er'^ + (-12 + QE + 2E'^ - E^y]Dl^ + aD^ + /3 = 0, (12) 

where a and /3 depend on r and E, naturally. For E"^ = 4 
this polynomial, Eq.(12), reduces to degree 2 instead of 3, 
and the multivalued Dec has indeed first order poles for 
E = 2 and E = —2, with residues 1 and r^, respectively. 
These agree with the exact residues, see Eq.(4). 

Finally on Fig. 6 we compare, for r=l.l, the exact D, 
as a function of E, with two kinds of estimates obtained 
with c = 1. The small-dashed line indicates the real root 
of Eq.(12), already found to show poles a.t E = ±2. The 
long-dashed line shows the real part of the other two roots 
of Eq.(12). Its complete lack of pole structure makes it 
much less valuable than the real root. 

In conclusion, this letter gives three results for the 
variational theory of propagators. The first one is the 
existence of a non trivial and soluble model for mean 
field approximations, which are notoriously non linear 
and demand a sorting of their (intricate) solutions. The 
second result is indeed a complete classification of the so- 
lution branches of the model, with their physical mean- 
ing, namely symmetry conservation or breaking. We can 
keep track of their interconnection or lack of connection. 
We can also assess the validity domain of the amplitude 
estimate given by each branch. Of special interest are the 
time reversal breaking solutions which we discovered, be- 
cause they can provide the exact amplitude and exist only 
if enough correlation is included in the trial states. When 
correlations are excluded, parity breaking solutions still 
exist and are absolutely necessary to account for several 
singularities of the exact amplitude. Fully symmetric so- 
lutions are not flexible enough, in the uncorrelated vari- 
ational space, to account for all the singularities of the 
exact, fully symmetric solution, which is unique and cor- 
related. The third result is then the necessity of a linear 
admixture of such solution branches, to remove any am- 
biguity of choices between them, patch their limited va- 
lidity domains, and restore the uniqueness of the physical 
solution. As a first step for the representation of physi- 
cal correlations, mean field solutions can thus be better 
analyzed along the lines described by this model. 
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FIG. 1. Contour plot of F, sec Eq.(6), for {s,t) = {u,v) 
when (r, E, c) = (1, 3, —0.25). Vicinity of parity breaking sad- 
dle point {s,t) = (-1.19,0.13) with saddle value Fpb = 0.46. 



FIG. 2. Contour plot of F for {r,E,c) = (1, 1,-0.25) and 
(s,w) = {t,v). Vicinity of the time reversal breaking saddle 
point (s, u) = (-0.73, -0.068) with saddle value Ftb = 1.58. 



FIG. 3. Various estimates of amplitude D as functions of 

the constraint c for (r,E) — (—1.2, 1.5). Full lines: full sym- 
metry. Small-dashed line: parity breaking. Long-dashed line: 
time reversal breaking. Exact Dex = 0.51 at Cex = —1.29. 



FIG. 4. Fully symmetric estimates as functions of c for 
{r,E) = (1,1.5). Any folding of one of the full lines onto 
itself means a transition between two and four real estimates. 



FIG. 5. Amplitudes as functions of _B for r = —1.2. Full 
line: Dex- Dashed line: full symmetry estimate, c = 1. 
Short-dashed: pb, c = 1. Long-dashed: tb, c = 0. Poles de- 
mand full symmetry near E = ±2 and ph or th near E = 0. 



FIG. 6. Full symmetry estimates as functions of E for 

(r, c) = (1.1,1). Short dashes: F real. Long dashes: ReF, 
complex F. Full line: Dex ■ Poles E = ±2 recovered by F real. 
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